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Abstract 

Background: Ixodes ricinus is a major vector for a range of microbial pathogens and the most prevalent and widely 
distributed tick species on the European continent, occurring in both natural and urban habitats. Nevertheless, little 
is known about the relative density of ticks in these two ecologically distinct habitats and the diversity of tick-borne 
pathogens that they carry. 

Methods: We compared densities of questing /. ricinus nymphs and adults in urban and natural habitats in Central 
and Northeastern Poland, assessed the prevalence and rate of co-infection with A. phagocytophilum, Rickettsia, 
Ehrlichia and 'Ca. Neoehrlichia spp.' in ticks, and compared the diversity of tick-borne pathogens using molecular 
assays (PCR). 

Results: Of the 1325 adults and nymphs, 6.2% were infected with at least one pathogen, with 4.4%, 1.7% and 
less than 0.5% being positive for the DNA of Rickettsia spp., A. phagocytophilum, Ehrlichia spp. and Ca. N. 
mikurensis, respectively. Although tick abundance was higher in natural habitats, the prevalence of the majority 
of pathogens was higher in urban forested areas. 

Conclusion: We conclude that: (i) zoonotic genetic variants of A. phagocytophilum are widely distributed in the 
Polish tick population, (ii) although the diversity of tick borne pathogens was higher in natural habitats, 
zoonotic species/strains were detected only in urban forests, (iii) and we provide the first description of Co. N. 
mikurensis infections in ticks in Poland. 

Keywords: Ixodes ricinus, Tick, Anaplasma, Rickettsia, Ehrlichia, Neoehrlichia, Natural habitats, Urban habitats, Tick 
density, Prevalence 



Background 

Pathogens that cause tick-borne diseases (TBDs), such 
as borreliosis or anaplasmosis, constitute a significant 
problem for the health of humans, companion animals 
and livestock worldwide. The majority of TBDs are clas- 
sified as emerging infectious diseases because of the 
growing awareness of their importance, the availability 
of better diagnostic reagents and the consequent rising 
statistics for incidence in the human population [1,2]. In 
Europe, Ixodes ricinus is the most prevalent and widely 
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distributed tick species and serves as the most important 
vector for several microbial pathogens [3-5]. However, 
tick abundance and the prevalence of transmissible 
pathogens in ticks varies locally depending on many 
factors, i. e. temperature, accessibility of suitable verte- 
brate hosts etc. [6-9]. The range of available hosts dif- 
fers depending on the habitat characteristics (i.e. urban 
versus natural) and is likely to also influence the spe- 
cies/strain variability of pathogens vectored by ticks 
[10]. In urban habitats, humans, pets (mainly dogs), 
synantropic rodents {Apodemus sylvaticus) and birds 
probably play a significant role as tick hosts and the 
sources of TBPs. In natural forest habitats in Europe 
mainly deer (Cervidae), wild rodents (Myodes glareolus, 
Apodemus flavicollis) and birds are considered to be the 
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most important hosts for /. ricinus and the major source 
of manyTBDs [11-14]. 

In Poland, studies on the occurrence and diversity of Rick- 
ettsia, Anaplasma, Ehrlichia and 'Candidatus Neoehrlichia' 
in I. ricinus have been few and the subject has been rela- 
tively neglected when compared with other TBPs [15,16]. 
The best recognized is the occurrence of A. phagocytophi- 
lum in /. ricinus ticks, but little is known about the existence 
of co-infections and the species/strain diversity of other 
Rickettsiaceae and Anaplasmataceae species. 

Tick-borne rickettsiosis is caused by intracellular bacteria 
belonging to the spotted fever group (SFG) of the Rickett- 
siaceae family. Since the early 1990s and the introduction 
and increasingly popular use of molecular methods (PCR, 
PCR-RFLP and sequencing), 17 new Rickettsia species have 
been added to the list of rickettsiae that are known to cause 
human diseases in different parts of the world [17]. Ticks 
are believed to act as vectors and reservoirs of the SFG 
group Rickettsia and interestingly, Rickettsia spp. are 
among the most common pathogens found in Dermacen- 
tor reticulatus and /. ricinus ticks in Europe [18]. However, 
reservoir hosts of Rickettsia are still not well defined and 
further investigations are required to elucidate the full 
range of potential host species. 

The family Anaplasmataceae are gram-negative, intra- 
cellular bacteria, and include the genera Anaplasma and 
Ehrlichia as well as bacteria of the new candidate species 
'Candidatus Neoehrlichia mikurensis' and 'Candidatus 
Neoehrlichia lotoris'. It is already well established that 
/. ricinus ticks constitute competent vectors for A. 
phagocytophilum, an agent of human granulocytic ana- 
plasmosis [19]. These bacteria are widespread over the 
territories of the USA, Europe and Asia and, apart from 
humans, can infect many mammalian species, i.e. dogs, 
cats, horses, sheep, rodents, birds and roe deer [14,20-22]. 
Molecular analyses have indicated that some strains/ gen- 
etic variants of A. phagocytophilum that are pathogenic 
for humans and domestic animals, circulate widely in na- 
ture in different hosts and display different vector tro- 
pisms and degrees of host pathogenicity [23-26] . 

Bacteria of a new candidate species 'Candidatus 
Neoehrlichia mikurensis' were discovered for the first 
time recently in small rodents and /. ovatus ticks in 
Japan [27]. It is believed that small rodents serve as 
reservoir hosts for Ca. N. mikurensis [28,29]. More- 
over, human infections with these bacteria have been 
described recently in North and Central Europe as well 
as in China [30-33]. The occurrence of Ca. N. mikur- 
ensis DNA in /. ricinus ticks has been confirmed already 
in western and northern parts of Europe [34-36] and in 
/. persulcatus ticks in eastern Russia [20]. 

Among the five described Ehrlichia species, three 
cause infections in dogs and humans (E. canis, E ewingii 
and E. chaffeensis), E. maris is a rodent pathogen and 



E. ruminantium is an agent of heartwater disease in 
domestic ruminants [20]. The first case of E. muris 
infection was described originally in a wild-caught 
specimen of the mouse Eothenomys kageus in Japan 
[37] and since then single cases of E. muris infection 
have been detected in rodents [38] and /. ricinus ticks 
[39] in Europe. In Poland only one report has been 
published thus far and that was for the occurrence of 
Ehrlichia sp. in common voles, Microtus arvalis [40]. 

The key hypothesis driving our research was that tick 
abundance and hence the prevalence and the diversity of 
the TBPs that they carry should differ markedly between 
natural and urban forests where ticks abound, because 
of the contrasting ecology with either little or consid- 
erable human intervention, respectively, and the avail- 
ability of different hosts for the ticks and hence 
reservoirs of varying competence for the pathogens. 
Our aims were therefore (1) to compare the densities 
of questing /. ricinus ticks (nymphs and adults) in for- 
ested urban and natural habitats; (2) to estimate the 
prevalence of A. phagocytophilum, Rickettsia, Ehrlichia 
and 'Ca. Neoehrlichia spp.' in /. ricinus ticks in these 
two different habitat types (urban and natural), and (3) 
to compare the species/strain diversity of tick-borne 
pathogens in two representative habitats in Central 
and Northeastern Poland. 

Methods 

Study sites 

Five study sites in Central and North Eastern Poland 
were investigated in 2011. The two sites selected as ex- 
amples of urban areas are located within the adminis- 
trative borders of Warsaw city, near the city center 
(<8 km). Natural areas are situated in protected areas 
of national or landscape parks in rural settings at some 
distance from major conurbations. 

Urban areas 

The Bielahski [5217'32"N, 20°57'36"E] and Kabacki Forests 
[52°6'58"N, 21°3'26"E] act as large city parks (130 ha 
and 903 ha, respectively), surrounded by dense urban 
and sub-urban settlements and are fragmented by paved 
roads and pedestrian/ cycle paths (Figure 1). The Kardynal 
Wyszyhski University is situated in Bielahski Forest. 
Access to urban forests is facilitated by main roads, the 
metro line and is within easy walking distance from 
nearby housing estates. These forests are affected by 
high anthropopression, being used daily as recreational 
areas by children, and intensively throughout the year 
by amateur cyclists, runners and country ski runners. 
The forests are also frequented by free-ranging feral 
dog populations and cats from the local housing es- 
tates. In Kabacki Forest, visitors are permitted to leave 
the paths to access facilities such as playgrounds for 
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Figure 1 Map of the natural and urban areas in Central, North-eastern and East Poland where ticks were collected. 



children, sports grounds, banks, approved camp fire/bar- 
beque sites. In these urban parks, tick collection sites were 
selected near to the main park entrance (<1 km) and main 
roads. The fauna of Kabacki and Bielahski forests is mainly 
limited to synantropic mammals such as rodents, urban 
foxes, hedgehogs, martens, badgers and a variety of spe- 
cies of birds. Roe deer, hare and wild boar may also occa- 
sionally venture into these locations. Elk migrating from 
the Kampinoski National Park may be occasionally en- 
countered in Bielahski Forest. 

Natural areas 

Three of the five study sites (Kampinoski National 
Park, KNP [5218'21"N, 20°36'32"E], Biatowieza Na- 
tional Park, BNP [52°46'20"N, 23°50'60"E], Mazury 
Landscape Park, MLP [53°48'25"N, 21°38'36"E]) are 



situated in protected areas of national or landscape 
parks and were considered as 'natural' areas, in Central 
(KNP) and Northeastern Poland (BNP, MLP) (Figure 1). 
Biatowieza National Park (62 500 ha) has a global 
reputation as a model natural habitat, representative of 
the ancestral European primeval deciduous forests 
(World Heritage List of UNESCO since 1972). It is 
known also for its populations of game species, includ- 
ing roe and red deer and wild boar, and in particular 
for its free-living population of European bison (Bison 
bonasus). Among other mammalian species regularly 
encountered in the forest are wolves, elks, lynx and 
beavers. 

Both Kampinoski National Park and Mazury Landscape 
Park comprise a mixture of natural and managed forests, 
some of the latter having been managed for over 120 years 
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in the case of MLP. Managed forests consist of pine trees 
{Pinus sylvestris) with the addition of deciduous tree 
species (birch, alder, elm and oak). MLP is well known 
for its combination of impressive forests and numerous 
(>1000) natural lakes. These forests are not fragmented, 
but are surrounded by agricultural areas, with minimal 
numbers of human settlements and are inhabited by a 
range of wild game and protected species indigenous to 
the region. As with BNP these include: roe and red 
deer, wild boar, elk, lynx, wolf, beaver and many other 
species of small and medium size mammals (i.e fox, 
badger, hare, marten) and birds. Human activity in the 
forests is very limited and seasonal; hunters, local people 
collecting berries or mushrooms, tourists (mainly cyclers), 
forestry workers. It is mandatory to keep dogs on the leash 
in the national parks in Poland, otherwise entry is prohib- 
ited by law and contravention is enforced by fines. 

As far as we are aware there are no quantitative census 
data through which anthropopression can be reliably 
compared between our urban and natural study sites, 
but from our personal experience of having visited our 
study sites regularly we estimate the number and density 
of visitors to be 1000-10000 higher in the urban parks 
compared to the natural forests. Similarly, the predator 
pressure from domestic and feral dogs and cats is un- 
doubtedly much higher in the two urban forests, as op- 
posed to the natural predator-prey relationships that 
exist in natural forests. 

Ixodes ricinus ticks 

Questing ticks, only nymphs and adults, were collected 
from vegetation by blanket dragging (white blanket 
1.5x0.75 m) in central and north-eastern regions of 
Poland in May, July and September 2011. Ticks were 
collected twice daily (about 9 a.m. and 4-5 p.m.), for a 
minimum of four consecutive days in each month. In 
each collection site, tick abundance was measured by 
conventional flagging in 300-500 m 2 plots and then ex- 
trapolated to a standard 100 m 2 . The species, develop- 
mental stage of individual ticks, as well as the sex of 
adults were recorded, and the nymphs from a single drag 
during flagging were pooled in groups of 10 prior to 
DNA extraction. In the case of positive PCR results for 
nymphs, it was assumed that at least one nymph was in- 
fected (minimal infection rate; MIR). Genomic DNA was 
extracted from homogenized tick specimens using the 
DNAeasy Blood and Tissue Kit (Qiagen, Crawley, UK ) 
and stored at -20°C. 

PCR analysis 

Detection and genotyping of Rickettsia, A. phagocytophi- 
lum, Ehrlichia spp. and Ca. Neoehrlichia spp. were per- 
formed by amplification and sequencing of different loci: 
the groESL heat shock operon and 16S rRNA gene for 



the Anaplasmataceae family and gltA for the Rickettsia- 
ceae family. The primers and thermal profiles used in 
this study have been described previously [19,41,42]. Re- 
actions were performed in a final volume of 20 ul and 
contained 0.33 mM dNTPs (Eurobio, Lille, France), 
2 mM MgCl 2 , lx PCR buffer, 1 U Taq polymerase (Fer- 
mentas), 1 uM of each primer and 5 ul of the extracted 
DNA sample. Negative controls were performed in the 
absence of template DNA. To check whether amplifiable 
DNA had been extracted, we used PCR reactions employ- 
ing tick-specific primers [43]. Only positive samples were 
chosen for further analysis. Amplicons were visualized 
with Midori Green stain (Nippon Genetics Europe 
GmbH) following electrophoresis in 2% agarose gels. 
Amplicons were purified using the Axygen Clean-up 
purification kit (Axygen, USA) and sequenced by a private 
company (Genomed S.A., Poland) in both directions. 

Phylogenetic analysis 

DNA sequence alignments and phylogenetic analysis 
were conducted using MEGA version 5.0 [44]. After test- 
ing the data for the best substitution model, phylogen- 
etic trees were obtained using Maximum Likelihood as 
the tree construction method and Hasegawa-Kishino- 
Yano parameter algorithm as a distance method. For 
comparison, sequences of Rickettsia, A. phagocytophilum, 
Ehrlichia and Ca. Neoehrlichia species/ strains obtained 
from GenBank (www.ncbi.nlm.nih.gov) were implemented 
in the sequence alignment. The stability of inferred phy- 
logenies was assessed by bootstrap analysis of 1000 ran- 
domly generated sample trees. 

Statistical analysis 

Data are expressed as arithmetic means ± 95% confi- 
dence limits per 100 m of blanket dragging. Tick abun- 
dance was analyzed using multifactorial general linear 
models (GLM) with normal error structures in SPSS 
v.20 for Windows. We fitted quality of habitat (HABI- 
TAT, 2 levels- urban [Bielanski + Kabacki Forests] and 
natural [KNP + BNP + MLP]), SITE (5 levels: KNP, BNP, 
MLP, Bielanski Forest, Kabacki Forest) and tick STAGE 
of development (3 levels: nymphs, females and males) as 
factors. Differences in Rickettsia, A. phagocytophilum, 
Ehrlichia and Ca. Neoehrlichia prevalence in ticks be- 
tween study sites were analyzed using the test in 
Instat or Fisher Exact Test when applicable. 

New nucleotide sequences 

New nucleotide sequences have been deposited in Gen- 
Bank with the accession numbers KF312351-KF321354 for 
16S rRNA of A. phagocytophilum, KF312355-KF321361 for 
groESL of A. phagocytophilum, KF312362 for groESL of 
E. maris and KF312363 for groESL of Ca. N. mikurensis. 
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Results 

Abundance of /. ricinus ticks 

Overall a total of 1589 ticks were collected and of these 
72% (1147) were nymphs and 28% were adults (241 
males and 201 females) (Table 1). The abundance of 
questing /. ricinus (nymphs and adults) varied signifi- 
cantly between study sites (GLM, main effect of SITE; 
F 4 ,i4o = 4.80, P = 0.001). The highest tick density per 
100 m 2 was recorded in Biatowieza National Park 
(Table 1). Average tick densities were similar in Bielahski 
Forest and MLP. The overall abundance of I. ricinus was 
lowest in KNP and Kabacki Forest (Table 1). The abun- 
dance of nymphs was highest in Bielahski Forest and 
BNP (GLM, main effect of SITE; F 4 , 140 = 3.90, P = 0.002, 
Table 1). Female and male abundances were relatively 
low (0.6-2.8 ticks per 100 m 2 ) and no significant differ- 
ences were observed between sampling areas. The abun- 
dance of /. ricinus ticks was almost twice as high in 
natural as in urban areas (GLM, main effect of HABI- 
TAT; F ljl40 = 5.80, P = 0.017) (Table 1). Because the dif- 
ference in the average tick densities per 100 m 2 between 
the two types of forest was significant, the sampling 
areas were divided into two groups with: (1) low tick 
abundance (5.7 [2.7-8.7]; Kabacki Forest and KNP), and 
(2) high tick abundance (14.5 [12.2-16.7]; Bielahski 
Forest, BNP, MLP) (GLM, main effect of HABITAT; 
Fi,i4o = 2L2, P< 0.001) (Table 1). No statistical differ- 
ences in tick abundance between study months were 
observed. 

Prevalence of pathogens in ticks 

Of the 1325 adults and nymphs, a total of 6.2% (« = 82) 
were infected with at least one pathogen. As expected, 
the prevalence of pathogens in ticks was significantly 
higher in adults (10.0%; 42/421) than the MIR in nymphs 
(4.4%; 40/904) Of 2 = 14.31, df=l, P< 0.0001). Prevalence 
of pathogens in ticks differed between urban and natural 
habitats (Table 2). Overall prevalence of TBPs was almost 
three times higher in urban compared with natural forests 
(X 2 = 22.97, df=l,P< 0.0001). 

Rickettsia spp. infections 

This was the predominant pathogen in our study, 4.4% 
(58/1325) of the tested ticks being recorded as positive 
for Rickettsia spp. The prevalence of infected ticks was 
about 2.5 times higher in urban (7.7%) than in natural 
sites (2.9%) (j = 7.77, df= 1, P = 0.0053). The DNA of 
Rickettsia spp. was detected in 3.7% nymphs (MIR) (33/ 
904) and in 5.9% adults (25/421) (NS). Prevalence of 
Rickettsia was significantly higher in low-tick-density 
forests (7.8% in KNP + Kabacki Forest) than in high-tick- 
density forests (3.1% in BNP + MLP + Bielahski Forest) 
(j 2 = 12.60, df= l,P< 0.001). 



Sequence analysis of the 770 bp fragment of the gltA 
gene of 40 isolates (20 of each habitat type, urban and 
natural) showed the presence of two different Rickettsia 
species. The most prevalent was R. helvetica (38/40, 
95%), closely related (99.9% homology) to the R. helve- 
tica isolates [GenBank: EU359285] originally obtained 
from Ixodes ticks in Switzerland and differed from them 
by only one nucleotide (G — > T) at position 713. Rickettsia 
helvetica isolates were found in all study sites. Two of 
40 isolates (5%) were identified as R. monacensis. These 
isolates showed 99.7% similarity (nucleotide substitu- 
tions at positions 476 [T -> C] and 713 [G->C]) to R. 
monacensis IrR/Munich strain [GenBank: DQ100163] 
derived from /. ricinus in Germany and were identical 
with R. monacensis isolated from human patients in Spain 
[GenBank: DQ517498] and Korea [GenBank: FJ009429]. 
The two R. monacensis isolates both originated from an 
urban forest (Kabacki Forest). 

Anaplasma phagocytophilum infections 

The overall prevalence of A. phagocytophilum in /. rici- 
nus was 1.7% (22/1325). Prevalence was significantly 
higher in urban than in natural sites (3% versus 1.1%) 
(Table 2, x = 4.98, df=l, P= 0.026). No positive ticks 
were found in Kabacki Forest and Mazury Landscape 
Park forest. The MIR in nymphs was significantly lower 
(0.6%; 5/904) than the prevalence in adults (4.0%; 17/ 
421) (x 2 = 18.26, df= 1, P < 0.0001), and prevalence did 
not differ statistically between ticks from low- and high- 
tick-density forests (Table 2). 

The 540 bp fragment of the 16S rRNA gene and the 
1200 bp fragment of the groESL heat shock operon were 
further analyzed in 18 isolates. The nucleotide identity/ 
similarity of the sequenced 16S rDNA fragments was 
very high (99.6-100%). Eight of 18 sequences from nat- 
ural areas (Table 3; KNP + BNP) were identical, repre- 
senting genetic variant I, which is considered to be 
zoonotic (Figure 2). A further eight isolates could be dis- 
tinguished on the basis of substitution at position 376 
(G — > A) in the variable region near the 5' end of the 16S 
rRNA gene, representing genetic variant II (Table 3). 
Variant III was represented by one isolate and differed 
by two nucleotides at positions 199 (G — » T) and 376 
(G — » A). Variant IV included one isolate that differed by 
two nucleotides at position 376 (G — > A) and 539 (G — ► T). 
Isolates belonging to variant I were found in ticks collected 
from natural areas. Isolates of variants II, III and IV were 
found only in ticks from one urban forest- Bielahski Forest 
(Table 3). 

The nucleotide sequences of isolates from groups I 
and II were identical to the sequences obtained from 
rodents, dogs, cats, roe deer, hedgehog and ticks 
(Table 3). Additionally, the 16S rRNA sequence of 
variant I was identical with the human pathogenic 



Table 1 Abundance of /. ricinus ticks in forested urban and natural regions of Poland 



Collection Quality 
site of habitat 



Abundance of Ixodes ricinus ticks (per 100 m [95% CL]) 



Total, all months 



May 



July 



September 





Nymphs 


Females 


Males 


Total 


Nymphs 


Females 


Males 


Total 


Nymphs 


Females 


Males 


Total 


Nymphs 


Females 


Males 


Total 


Bielariski urban 


8.3 


0.9 


1.2 


10.5 


17.5 


4.0 


0.5 


22.0 


9.5 


2.0 


1.5 


13.0 


11.8 


2.3 


1.1 


15.2* 


8.8 


Forest 


[5,3-11,3] 


[0,1-1,8] 


[0,2-2,3] 


[6,8-14,1] 


[5,4-29,6] 


[0,4-7,6] 


[0,0-4,9] 


[7,2-36,8] 


[0,0-21,6] 


[0,0-5,6] 


[0,0-5,9] 


[0,0-27,8] 


[6.0-1 7.6] 


[0.6-4.0] 


[0.0-3.2] 


[8.1-22.2] 


[5.9-11.6] 


Kabacki 


3.7 


02 


0.3 


4.2 


6.0 


1.5 


1.0 


8.5 


95 


0.0 


15 


11.0 


6.4 


0.6 


0.9 


7.9 




Forest 


[0,0-7,5] 


[0,0-1,3] 


[0,0-1,6] 


[0,0-8,8] 


[0,0-18,1] 


[0,0-5,1] 


[0,0-5,4] 


[0,0-23,3] 


[0,0-21,6] 


[0,0-3,6] 


[0,0-5,9] 


[0,0-25,8] 


[0.6-12.2] 


[0.0-2.3] 


[0.0-3.0] 


[0.8-15.1] 




BNP natural 


115 


3.4 


3.9 


18.8 


21.3 


2.3 


2.7 


26.3 


8.2 


2.2 


1./ 


12.1 


13.7 


2.6 


2.8 


19.1* 


13/1 




[7,6-15,5] 


[2,2-4,5] 


[2,5-5,4] 


[14,0-23,6] 


[14,3-28,3] 


[0,3-4,4] 


[0,2-5,2] 


[17,8-34,9] 


[2,5-13,9] 


[0,5-3,9] 


[0,0-3,7] 


[5,1-19,1] 


[10.4-17.0] 


[1.7-3.6] 


[1 .6-3.9] 


[15.1-23.1] 


[10.9-16.0] 


KNP 


1,7 


2.1 


2.9 


6,8 


2.3 


0.8 


1.3 


4.5 


6.0 


0.3 


1.0 


7.3 


3.4 


1.1 


1.8 


6.2 






[0,0-6,2] 


[0,8-3,4] 


[1,3-4,5] 


[1,4-12,2] 


[0,0-9,3] 


[0,0-2,9] 


[0,0-3,8] 


[0,0-13,0] 


[0,0-14,6] 


[0,0-2,8] 


[0,0-4,1] 


[0,0-1 7,7] 


[0.0-7.3] 


[0.0-2.2] 


[0.3-3.2] 


[1.3-11.0] 




MLP 


13.3 


0.3 


0.5 


14,1 


3.0 


0.0 


5.0 


8.0 


10.0 


2.0 


2.4 


14.4 


8.8 


0.8 


2.6 


12.2* 






[8,6-18,1] 


[0,0-1,7] 


[0,0-2,2] 


[8,3-19,9] 


[0,0-20,1] 


[0,0-5,0] 


[0,0-11,2] 


[0,0-28,9] 


[2,3-1 7,7] 


[0,0-4,3] 


[0,0-5,2] 


[5,0-23,8] 


[2.3-15.2] 


[0.0-2.7] 


[0.3-4,9] 


[4.2-20.0] 




Total, all sites 


/./ 


1.4 


1.7 


10.9 


10.0 


1.7 


2.1 


13.9 


8.6 


1.3 


1.6 


11.6 


8.8 


15 


1.8 


12.1 [9.3-14.9] 




[5.9-9.5] 


[0.9-1.9] 


[1.1-2.4] 


[8.7-13.1] 


[4.8-15.2] 


[0.2-3.3] 


[0.2-4.0] 


[7.5-20.3] 


[4,4-12.9] 


[0.0-2.6] 


[0.1-3.2] 


[6.3-16.8] 


[6.5-11.1] 


[0.8-2.2] 


[1 .0-2.7] 







CL- 95% confidence limits, 
"classified as high-tick-density forests. 
N- nymphs, F- females, M-males. 
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Table 2 The prevalence of tick-borne pathogens in /. ricinus collected in forested urban and natural areas of Poland 



Site Quality Overall TBD 

of habitat prevalence% 



Prevalence% (no. infected/no. tested) 



Rickettsia sp. 



A. phagocytophilum Ca. Neoehrlichia sp. 



Ehrlichia sp. 



Bielariski Forest urban 


1 1 .6 (42/405) 


4.5 (9/201) 


7.7 (31/405) 


6.0 (12/201) 


3.0 (12/405) 


0.0 (0/201) 


0.2 (1/405) 


0.0 (0/201) 


0.7 (3/405) 


Kabacki Forest 




1 0.8 (22/204) 




0.0 (0/204) 




0.5 (1/204) 




1 .5 (3/204) 




Biatowieza natura 
National Park 


4.4 (40/920) 


1.6 (10/625) 


2.9 (27/920) 


1 .3 (8/625) 


1.1(10/920) 


0.0 (0/625) 


0.2 (2/920) 


0.0 (0/625) 


0.1 (1/920) 


Kampinos 
National Park 




3.8 (6/1 56) 




1 .3 (2/1 56) 




0.0 (0/1 56) 




0.0 (0/156) 




Mazury Lake 
District* 




7.9 (11/139) 




0.0 (0/139) 




1 .4 (2/1 39) 




0.7 (1/139) 





"classified as high-tick-density forests. 

strains isolated in Slovenia, Italy and the USA (Table 3). 
The nucleotide sequences of isolates from groups III 
and IV did not show 100% identity to the known se- 
quences of 16S rRNA of A. phagocytophilum deposited 
in GenBank. 

The partial groESL fragments (1200 bp) were also se- 
quenced for all positive samples (n = 18) since the 16S 
rRNA gene is too conserved for analysis of genetic het- 
erogeneity. The level of homology between isolates was 
also high (98.0-100%). Sequence analysis allowed identi- 
fication of seven different genetic variants (I-VII) that 
differed by 2 to 13 nucleotides (Table 4). Genetic vari- 
ants I and II included 10 isolates from ticks collected 
only in Bielahski Forest, an urban site (Table 4). Variants 
III-VII were composed of 1 or 2 sequences found only in 
ticks from natural forested areas (BNP, KNP). Analysis 



of these groESL fragments allowed the differentiation of 
4 variants from BNP, all grouped in variant I based on 
the 16S rRNA gene fragment (Table 3). 

Scrutiny of the phylogenetic tree, based on the partial 
groESL operon sequences, showed that two isolates from 
variant VII were identical to other European isolates, 
mainly from /. ricinus ticks and roe deer (Figure 2, 
Table 4). Isolates belonging to variant Till clustered 
with A. phagocytophilum pathogenic for human and 
domestic animals in Europe, as well as in North America, 
and differed from variant VII by 13 nucleotides (Table 4). 
Thus, these 6 variants (I- VI) are considered to represent 
zoonotic genotypes. The nucleotide sequences of isolates 
from groups IV- VI did not show 100% identity to the 
known sequences of groESL operon sequences of A. pha- 
gocytophilum deposited in GenBank. 



Table 3 Polymorphism in the fragment of 16S rRNA gene in A. phagocytophilum isolates from ticks and human 
pathogenic strains (sequences published in GenBank) 



Strain/genetic 
variant 



Nucleotide 
positions 5' — » 3' a 

199 328 376 539 



Quality of habitat/ 
country 



No. of isolates 
(reference no. of 
samples)/site 
of study 



GenBank Host/vector found in other studies 
acc. (only 100% homology; nucleotide 

number sequences deposited in GenBank) 



This study 



natural/ Poland 



2 (6,153)/KNP 

6 (39, 54, 73, 85, 5 
283)/BNP 



KF31 2352 bank vole (Myodes glareolus, KC740432); 

dog (KF985242); elk (Alces alces, KC800983); 
roe deer {Capreolus capreolus, JQ965530); 
cotton rat {Sigmodon hispidus, JQ063025); 
rat (Rattus norvegicus, KC470064); reindeer 
{Rangifer tarandus, JX841254); hedgehog 
{Erinaceus europaeus, JN571 163); /. ricinus 
(KF481930), /. trianguliceps (KF481934) 



II 


G 


A 


A 


G 


urban/Poland 


8 (29, 37, 38, 47, 62, 67, 
71, 73)/Bielahski Forest 


KF3 12353 


dog (EF668225); cat (HM1 38366); /. ricinus 
(JX909354); D. reticulata (KF381413) 


III 


T 


A 


A 


G 




1 (57)/Bielahski Forest 


KF3 12351 


n. m. b 


IV 


G 


A 


A 


T 




1 (64)/Bielahski Forest 


KF3 12354 


n. m. b 


A. phagocytophilum 
strains pathogenic 
for human 


G 
G 
G 


A 
A 
A 


G 
G 
G 


G 
G 


Poland 
Slovenia 
Italy 




KF111754 
GU236658 
DQ029028 






G 


A 


G 


G 


USA 




GU236664 






G 


A 


G 


G 


USA 




NRJJ74113 




a The number corresponds to the positions of nucleotide substitutions relative to the sequence of the 
(NR_074113). Base substitutions are shown in bold. 


complete 16S rRNA gene of A. phagocytophilum strain HZ 



No match has been found in GenBank. 
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Variant II 

EU860089 sheep France 
Variant V 
Variant III 

■ HM057224 red deer Spain 
Variant VI 

HM057232 Ixodes ricinus Spain 
AF202895 Ixodes ricinus Switzerland 
AF548385 sheep Norway 
Variant IV 
Variant I 

EU381151 A p dog Slo\enia 
EU381 152 A p Ixodes ricinus Slovenia 
HQ629910 Ixodes ricinus Estonia 
AF033101 human Slo\enia 
EU381150 dog Slo\*nia 
I JF494839 human USA 
I U96728 human USA 
AF172159 human USA 
AF172163 human (USA) 
JF494833 dog USA 
AY281830 Ixodes ricinus Germany 
AY220469 roe deer Austria 
AF478556 roe deer Sloysnia 
AY281811 Ixodes ricinus Germany 
HQ629905 Ixodes ricinus Estonia 
AF478560 roe deer Slo\enia 
AF383225 roe deer USA 
Variant VII 

HQ629902 A p Ixodes ricinus Belarus 
AY220468 roe deer Czech Republic 
AY281817 Ixodes ricinus Germany 
AY281850 A p Ixodes ricinus Germany 
DQ779568 roe deer Poland 
EU246960 Ixodes ricinus Slovenia 
JN005746 roe deer Poland 
AY044161 A platys 

our Ehrlichia 

GU358690| E. muris isolate Sv-40 
GU358687 E. muris isolate Nov-lp300 
GU358689| E. muris isolate Kh-1550 
JQ085941 E. chaffeensis 
ECU96731 E. canis 



92 

96 

I AF 



Anaplasma phagocytophilum 



- AF195273 E. ewingii 
- DQ647009 E.a ruminantium 

EF633745 Ca. N. lotoris racoon rocyon lotor USA 

EU810407 Ca. N. mikurensis tick Ixodes ricinus Germany 
JQ669946 Ca. N. mikurensis Ixodes ricinus Italy 
EU810406 Ca. N. mikurensis human Homo sapiens Germany 
EU432375 Ca. N. mikurensis dog Canis familiaris Germany 
our Ca N mikurensis 

JX392803 Ca. N. mikurensis Myodes glareolus France 
FJ966363 Ca.N. mikurensis tick Ixodes persulcatus Russia 
JQ359062 Ca. N. mikurensis human Homo sapiens China 
JQ359070 Ca. N. mikurensis Rattus norvegicus China 
JQ359077 Ca. N. mikurensis Ixodes persulcatus China 
KC108716 Ca. N. mikurensis Apodemus agrarius China 
KC10S71 7. 1 1 Ca. N. mikurensis Apodemus peninsulae China 



Figure 2 Phylogenetic tree of the A. phagocytophilum, Ehrlichia i Ca. Neoehrlichia isolates studied in the current work and chosen 
isolates from GenBank based on the fragment of the groESL heat shock operon. Numbers at the nodes of the tree indicate bootstrap 
values (1000 replicates). 



75 FT 

> "Tl 

Table 4 Polymorphism in the fragment of the groESL heat shock operon in A. phagocytophilum isolates from ticks and human pathogenic strains (sequences published | Jj- 

in GenBank) 



Strain/genetic Nucleotide positions 5' — > 3 ,a 

variant 





167 


398 


401 


431 


476 


486 


497 


51 1 


561 


776 


821 


836 


837 


840 


845 


857 


872 


1 


G 


A 


A 


G 


G 


G 


T 


G 


A 


G 


G 


A 


c 


T 


A 


G 


c 


II 


G 


A 


A 


G 


c 


G 


c 


c 


A 


c 


c 


A 


c 


T 


A 


c 


c 


III 


G 


A 


A 


A 


c 


G 


T 


c 


A 


c 


c 


A 


c 


T 


A 


c 


c 


IV 


G 


A 


A 


G 


A 


A 


T 


G 


A 


c 


c 


A 


c 


T 


A 


c 


c 


V 


G 


A 


A 


G 


C 


G 


C 


C 


A 


c 


c 


A 


c 


T 


A 


c 


c 


VI 


A 


A 


A 


G 


c 


G 


C 


c 


A 


c 


c 


A 


c 


T 


A 


c 


c 


VII 


G 


G 


A 


G 


c 


G 


T 


c 


A 


c 


A 


G 


G 


G 


T 


T 


T 


im 


G 


A 


G 


G 


c 


G 


T 


c 


G 


c 


c 


A 


A 


T 


A 


c 


C 


liC 


G 


A 


A 


G 


c 


G 


T 


c 


G 


c 


c 


A 


A 


T 


A 


c 


c 




A 


A 


A 


G 


c 


G 


T 


c 


G 


T 


c 


A 


A 


T 


A 


c 


c 




G 


A 


G 


G 


c 


G 


T 


c 


G 


c 


c 


A 


A 


T 


A 


c 


c 



I — 

-S to 

V 

to * 

QJ «-+ 

QJ Q 

II 

S ^ 

g go 

n fV 

o H 

1 1 

n «< 

O ISJ 

13 O 

H ? 

NJ 



This study 



for human 



Strain/genetic 










Nucleotide positions 5' 


>3' a 








Quality of 


No. of isolates (reference 


GenBank acc. 


Host/vector found in other studies 


variant 


890 


923 


956 


965 


980 


992 


1049 


1085 


1088 


1112 


1229 


1244 


habitat/country 


no. of samples)/site of 
study 


number 


(only 100% homology; nucleotide 
sequences deposited in GenBank) 


This study 


G 


C 


G 


G 


T 


T 


C 


A 


C 


C 


A 


T 


urban/Poland 


9 (29, 37, 38, 47, 57, 62, 67, 
71 , 73)/ Bieiariski Forest 


KF3 12357 


dog (EU381 151); /. nanus (EU381 152) 




T 


C 


T 


G 


T 


T 


T 


A 


c 


c 


G 


C 




1 (64)/ Bieiariski Forest 


KF3 12360 


red deer (Cervus elaphus, HM057225); I. 
ricinus (KF383239) 




T 


c 


T 


G 


T 


T 


T 


A 


c 


c 


G 


C 


natural/Poland 


2(85, 90J/BNP 


KF3 12358 


/. nanus (AY281 831) 




G 


c 


C 


G 


T 


T 


C 


A 


c 


c 


A 


T 




2(39, 54J/BNP 


KF3 12355 


n. m. b 




G 


c 


T 


G 


T 


T 


T 


A 


c 


c 


A 


T 




K73VBNP 


KF3 12359 


n. m b 




G 


c 


T 


G 


T 


T 


C 


A 


c 


c 


G 


C 




1(283)/BNP 


KF312361 


n.m b 




G 


G 


C 


A 


A 


C 


T 


G 


c 


T 


G 


C 




2(6, 153)/KNP 


KF3 12356 


roe deer (Capreolus capreolus, 
AF478560); /. ricinus (EU552914) 


A phagoq/tophilum 


G 


c 


C 


A 


T 


T 


C 


G 


A 


T 


G 


C 


USA 




U96728 




strains pathogenic 
for human 


G 


c 


c 


A 


T 


T 


C 


G 


A 


T 


G 


c 


USA 




AF172159 






G 


c 


c 


G 


T 


T 


c 


A 


c 


C 


A 


T 


Slovenia 




AF033101 






G 


c 


c 


A 


T 


T 


c 


G 


A 


T 


G 


c 


USA 




JF494S39 





^The number corresponds to the positions of nucleotide substitutions relative to the sequence of the groESL heat shock operon of the human pathogenic A phagocytophilum strain (U96728). Base substitutions are shown in bold. 
b No match has been found in GenBank. 
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Ehrlichia spp. and Candidatus Neoehrlichia spp. infections 

The DNA of Ehrlichia spp. and that of Ca. Neoehrlichia 
spp. was found in 0.3% (4/1325; only in nymph pools) 
and 0.2% (3/1325; two males and 1 nymph pool) of the 
sampled ticks, respectively. Three ticks positive for Ehr- 
lichia DNA were collected in an urban area (Kabacki 
Forest) and one in MLP (natural site) (Table 2). The 
presence of Ca. Neoehrlichia spp. was confirmed in 
nymphs from Kabacki Forest (urban area) and in 2 males 
from MLP (natural site) (Table 2). 

Molecular analysis of the 540 bp fragment of the 16S 
rRNA gene of Ehrlichia (n = 4) and Ca. Neoehrlichia (n = 3) 
isolates showed that all sequences were identical and differ- 
entiation between these genera was possible only on the 
basis of the partial groESL fragment (1200 bp). The four Ehr- 
lichia isolates were identical [GenBank: KF3 12362] and 
showed 82.3% groESl sequence homology with Ca. Neoehr- 
lichia and 100% sequence homology with E. muris originally 
obtained from rodents [GenBank: GU358690] and I. persul- 
catus [GenBank: GU358686] in Russia. All Ca. Neoehrlichia 
isolates were also identical [GenBank: KF312363] and 
showed 100% sequence homology with Ca. N. mikur- 
ensis from /. ricinus ticks [GenBank: EU810407] and 
humans [GenBank: EU810406] in Germany. Our iso- 
lates clustered also with other Ca. N. mikurensis patho- 
genic for humans [GenBank: HM045824] and dogs 
[GenBank: EU432375] in Europe (Figure 2). 

Co-infections in /. ricinus ticks 

4.8% (2/42) of adult ticks identified as carrying pathogens 
yielded positive results for two pathogens and none were 
infected with three or more. Rickettsia helvetica and A. 
phagocytophilum were involved in both co-infections. 
Both ticks co-infected with these pathogens were col- 
lected in an urban area (Bielahski Forest). Analysis of 
co-infections in nymphs was not justified because DNA 
was extracted from pooled samples of ticks. 

Discussion 

The results of our study, focusing on a comparison of 
the abundance of /. ricinus ticks and prevalence of TBPs 
in ticks sampled in forested urban and natural habitats, 
revealed significant and interesting differences between 
these ecologically contrasting sites. Although tick abun- 
dance was higher in natural habitats, the prevalence of 
the majority of TBPs (A phagocytophilum, Rickettsia 
spp., E. muris) was higher in urban areas. Interestingly, 
while the diversity of A. phagocytophilum was higher in 
natural habitats, only the zoonotic species/ strains (i.e. R. 
monacensis, A. phagocythophilum variants I-II \groESL 
operon sequences]) were detected in urban forests. This 
finding likely reflects the availability of different reser- 
voir hosts and points to a higher risk of TBDs for humans 
venturing into urban forests compared with forests located 



in natural sites. We also report here for the first time the 
detection of Candidatus Neoehrlichia mikurensis in quest- 
ing Ixodes ricinus ticks in Poland. 

As in many previous studies on ticks [45,46], in our 
study also the abundance of questing /. ricinus varied 
significantly between study sites. The highest density 
was recorded in natural forested habitats, i.e. in the 
Bialowieza National Park and the Mazury Landscape 
Park both of which correspond to the 'medium tick 
abundance' category of forests (11-40 ticks per 100 m 
[47]). A similar high abundance was observed in Bielahski 
Forest in Warsaw, which is a part of an ancient deciduous 
primeval forest and the soil humidity and plant cover- 
age are very similar to those encountered in Bialowieza 
Primeval Forest (BNP). Deciduous humid forests are 
known to constitute preferable habitats for /. ricinus 
ticks [48,49]. The overall abundance of /. ricinus was 
lowest in Kampinos National Park and Kabacki Forest, 
where pine Pinus silvestris are the main tree species 
because of the poorer sandy soil conditions that are 
less suitable for deciduous tree species (mainly low- 
humid sandy areas; 'low tick abundance', 3-10 ticks per 
100 m 2 [47]). 

Despite lower tick abundance, the prevalence of three 
of four pathogens detected in this study was at least 
twice as high in ticks from urban forests, a finding cor- 
responding to the so-called 'dilution effect' of Estrada- 
Pena et al. [50]. The reasons underlying this difference 
and its consequences are of particular relevance for public 
health concerns. For an explanation of this phenomenon 
we looked at the genetic diversity of the most common 
pathogens. Although the diversity of A. phagocytophilum 
was higher in natural habitats, only zoonotic species/ 
strains (i.e. R. monacensis, A. phagocythophilum variants 
I-II \groESL operon sequences]) were detected in urban 
forests. This finding likely reflects the availability of differ- 
ent reservoir hosts for the pathogen and is consistent with 
the idea that in urban areas humans, companion animals 
(mainly dogs) as well as synanthropic rodents, foxes and 
hedgehogs constitute suitable hosts for ticks and act as 
reservoir hosts for the parasites. Rickettsia monacensis, 
which is pathogenic for humans, was found only in ticks 
from Kabacki Forest, supporting our hypothesis and 
pointing to a new health risk in this recreational area of 
Poland's capital city, Warsaw. The greatest genetic vari- 
ability of A. phagocytophilum (5 variants) was found in 
ticks from Bialowieza and Kampinos National Parks (in- 
cluding also the non-zoonotic variant VII \groESL operon 
sequences]) probably reflecting the wide range of natural 
I. ricinus hosts that live in these forests and can act as res- 
ervoir hosts for A. phagocytophilum (i.e. Cervidae, Carni- 
vores). In this context it is pertinent also that in our 
earlier studies on roe deer, the dominance of non- 
zoonotic over zoonotic A. phagocytophilum strains and 
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other TBPs (i.e. Babesia, Bartonella) was noted [14]. Thus, 
although tick density and TBP diversity in ticks in urban 
areas appeared to be lower than in the natural habitats 
and to depend on a narrower tick host range, the risk for 
humans may be higher, as humans 'share' many TBPs with 
dogs, livestock and rodents, as a result of many years of 
co-existence. 

Previous molecular studies have shown relatively high 
prevalence of tick-borne pathogens in /. ricinus from dif- 
ferent parts of Europe. The most common infection in the 
current study was Rickettsia spp. with a prevalence of 
4.4%, and this is comparable to prevalence rates reported 
from other countries in Europe: 5% in Luxembourg [51], 
14-17% in Germany [6], 1.4% in France [52], 11% in 
Slovakia [38] and 9% in Poland [53]. In contrast to the re- 
sults obtained by Overzier et al. [10] and Venclikova et al. 
[54], a significantly higher prevalence of Rickettsia was 
found in urban than in natural areas (7.7 vs. 2.9%). On the 
basis of the gltA gene we have identified two different 
Rickettsia species. Isolation of R. helvetica in our study 
confirmed the high prevalence and widespread nature of 
this species in /. ricinus in Europe. Since 1999, several 
R. helvetica seropositive patients with relatively mild, 
self-limited illnesses associated with flu-like symptoms 
have been reported, so the pathogenicity of R. helvetica. 
needs further investigation. The nucleotide sequences 
of the gltA fragment gene showed that our R. monacensis 
isolates were identical with the R. monacensis strain Rp- 
Spl that is known to be pathogenic for humans. R. mona- 
censis is an etiological agent of human rickettsiosis in 
Spain [55] and the DNA of these bacteria has been de- 
tected in an asymptomatic patient in Croatia [56]. The 
presence of these bacteria in /. ricinus ticks has also 
been confirmed in other European countries including 
Germany and Romania [57,58]. In Poland, R. monacen- 
sis has been detected recently in Northwestern Poland 
(in Szczecin, near Poland's western border with Germany), 
where this species has never been recorded previously in 
ticks [53]. 

A. phagocytophilum has been detected in /. ricinus ticks 
throughout Europe and the prevalence rate has been re- 
ported to vary from 1 to 5% in Poland, Switzerland, Russia 
and Belarus [59-61], however, locally it may exceed 24% 
(Italy [62]). As in our results with Rickettsia and the data 
of Venclikova et al. [54], we found a significantly higher 
prevalence of A. phagocytophilum in urban compared with 
natural areas (3.0 vs. 1.1%). In agreement with data pub- 
lished previously [26], sequence analysis of the 1200-bp 
fragment of groESL in the current work revealed two dis- 
tinct genetic lineages of A. phagocytophilum: (1) genetic 
variants detected in humans, ticks, dogs, horses, sheep 
and red deer from Europe and USA which are believed to 
be pathogenic; and (2) genetic variants isolated from ticks 
and roe deer in Europe that are probably non-zoonotic 



strains [61,63]. In our study, molecular analysis of the 16S 
rRNA fragment gene showed that the nucleotide se- 
quences of eight A. phagocytophilum isolates belonging to 
genetic variant I were identical with A. phagocytophilum 
strains known to be pathogenic for humans (Table 3). 
However, on the basis of molecular and phylogenetic ana- 
lysis of the groESL fragment, two isolates from variant I 
(no. 6 and 153), were assigned to the non-zoonotic variant 
VII. Therefore, our results emphasize that analysis of the 
nucleotide sequences of several genes is necessary, indeed 
essential, for a detailed delineation of different A. phagocy- 
tophilum strains [20]. 

This study has reported for the first time the presence 
of 'Ca. N. mikurensis' in questing /. ricinus from Poland. 
The prevalence rates of 'Ca. N. mikurensis' and E. muris 
in our study were low and differed from those in re- 
ports from Slovakia (E. muris 3% [39]), Switzerland 
(Ca. N. mikurensis 6% [64]) and Germany (Ca. N. 
mikurensis 6% [35]). In our study, the 540-bp fragment 
of the 16S rRNA gene was highly conserved and differ- 
entiation of the two genera was possible only on the basis 
of the groESL operon fragment. Although E. muris is 
mainly a murine pathogen, 'Ca. N. mikurensis' is now con- 
sidered to be an emerging tick-borne pathogen of veterin- 
ary and medical importance [65]. Isolates from this study 
were identical with 'Ca. N. mikurensis' detected in blood 
samples of patients with severe febrile illnesses and with 
the strain isolated from /. ricinus from Germany [32]. Our 
results confirm the presence of 'Ca. N. mikurensis' in 
questing /. ricinus from Poland and point to a new risk of 
infection with a recently discovered TBP for animals and 
humans. 

Interestingly, both ticks in which two genera of TBPs 
were detected (R. helvetica/ A. phagocytophilum) were 
collected in Bielahski Forests from the area of Warsaw. 
Thus, the risk of acquiring double infection from the in- 
fected ticks appears to be higher in urban areas. How- 
ever, because the number of detected co-infections was 
so low, on the basis of the current results it is not pos- 
sible to conclude whether these concurrently infected 
ticks can be considered as constituting an additional 
health risk for people visiting the city forests. 

Conclusions 

There were significant differences in tick abundance be- 
tween natural and urban forested habitats. Although tick 
abundance was higher in natural habitats, the prevalence 
of the majority of TBPs was higher in urban areas. This 
finding, together with the detection of zoonotic species/ 
strains of Rickettsiaceae and Anaplasmataceae indicates 
a high risk of TBDs in city forests. The new pathogen- 
Candidatus Neoehrlichia mikurensis- was detected for 
the first time in questing /. ricinus ticks in Poland, also 
in urban areas. 
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